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ABSTRACT 

Located within 10 degrees of the North Pole, northern Ellesmere Island offers continuous darkness in the 
winter months. This capability can greatly enhance the detection efficiency of planetary transit surveys and 
other time domain astronomy programs. We deployed two wide-field cameras at 80°N, near Eureka, Nunavut, 
for a 152-hour observing campaign in February 2012. The 16-megapixel-camera systems were based on com- 
mercial f/1.2 lenses with 70mm and 42mm apertures, and they continuously imaged 504 and 1,295 square 
degrees respectively. In total, the cameras took over 44,000 images and produced better-than- 1 % precision 
light curves for approximately 10,000 stars. We describe a new high-speed astro me trie and photometric data 
reduction pipeline designed for the systems, test several methods for the precision flat-fielding of images from 
very-wide-angle cameras, and evaluate the cameras' image qualities. We achieved a scintillation-limited pho- 
tometric precision of 1-2% in each 10s exposure. Binning the short exposures into 10-min chunks provided 
a photometric stability of 2-3 millimagnitudes, sufficient for the detection of transiting exoplanets around the 
bright stars targeted by our survey. We estimate that the cameras, when operated over the full arctic winter, will 
be capable of discovering several transiting exoplanets around bright (m v < 9.5) stars. 



1. INTRODUCTION 

The diurnal cycle is a significant impediment to ground- 
based transit surveys, limiting their detection efficie ncy for 
plane ts with orbital periods of tens of days or more (Brown 
2003 ). Consequently, exoplanet transit surveys rely on long- 
term observing campaigns and/or the use of multiple observ- 
ing sites to reach an acceptable probability of detecting mul- 
tiple transits. In contrast, circumpolar locations provide con- 
tinuous darkness in the winter months, potentially allowing 
single-instrument ground-based surveys to reach longer pe- 
riod planets with high detection efficiencies. 

The detec tions of tra nsiting pl anets around very bright stars 
(e.g. [Charbonneau et al. 2000; Henry et al. 2000; Sato et al. 



2005| [Bouchy et al. 2005; Christian et al. 2006; Bako s etal 
20071 |Pal et al.||2010||Winn et al.||20111|Howell et al.|2012| >, 
and the relative ease of characterization of those planets (e.g. 



Agol et al. 2010; Collier Cam eron et al.|2010||Demory et al. 
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von Braun et al.||201 1| |Stevenson et al.||201 1 Majeau 



2012| i, have demonstrated that very-small-aperture tele 



scopes can discover exciting exoplanets. Here, we describe 
two circumpolar very-wide-field, small-aperture cameras de- 
signed to search for transiting exoplanets around 5,000 stars 
brighter than my =9. 5. 

The advantages of circumpolar astronomy have already mo- 
tivated several site characterization and observational projects 
based on the high Antarctic glacial plateau. The Gattini cam- 
eras have measured the optical sky brightness, cloud cover 
and auroral emission at both Dome A and Dome C ( |Moore| 
et al.|2006| 12008} [20T0] >, while small, robotically controlled, 
optical and infrared telescopes are in the testing and operation 
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Fig. 1. — The High Canadian Arctic location of the test system (80°00'N, 
86°20'W) at the Ridge Lab near Eureka, Nunavut. The 2012 position of the 
North Magnetic Pole is also marked. 



phases fldi Rico et al.|2010[|Daban et al.|2010] >. The CSTAR 
project (Wang e t al.|201 1 1 demonstrated long-term photome- 
try on 10,000 stars in a 23 deg 2 region centered on the South 
Celestial Pole. 

In this paper we present the first astronomical survey con- 
ducted from Ellesmere Island, located in the High Canadian 
Arctic at 80°N in the territory of Nunavut (Figure [T]). Robotic 
site-testing systems deployed at nearby coastal mountain sites 
found a high clear-sky fraction, low median wind speed, and 
the possibility of excellent seeing (Steinbring et al. 20101. 
More intensive site testing at a lower-elevation manned sta- 
tion, the Polar Environment Atmospheric Research Labora- 
tory (PEARL), showed photometric conditions for approxi- 
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mately half of the continuously dark polar night, and continu- 
ous periods of hundreds of hours of clear, dark skies. Further- 
more, sufficiently clear conditions for spectroscopy or differ- 
ential photometry occur for up to 86% of the winter ( [Stein- 
bring et al.|2012| l. Seeing measurements at the same location 
show the presence of a weak turbulence layer near the ground 
( Hickson et al.|20f0| , potentially providing excellent observ- 
ing conditions. 

The excellent conditions at the Ridge Lab site and their im- 
plications for exoplanet detection efficiencies, along with the 
existing infrastructure and relative ease of accessibility, en- 
couraged us to design astronomical survey instruments for the 
site. In this paper we describe the design, operation and first 
results from two very-wide-field exoplanet-search cameras, 
the AWCam (Arctic Wide-Field Camera) prototypes, which 
we deployed to the Ridge Lab in February 2012. Although 
lasting only a week, and so within the regime for which a 
mid-latitude observatory could achieve similar efficiencies, 
this initial campaign was intended as a test of the survey capa- 
bilities achievable under real observing conditions at the site, 
as a prelude to longer operations in later years. 

The AWCam cameras continuously observed a several- 
hundred-square-degree field around the North Celestial Pole, 
using short exposures to avoid the need for a tracking mount. 
Like several lower-latitu de transi ting exoplanet s urveys such 
as SuperWASP 



( [Pollacco et al. 
2004), and KELT (Pepper et al. 2007 



2006), 



et al. 
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are based on commercial camera lenses. The cameras are 
primarily designed to search for transiting planets around 
bright (m v < 10) stars, but the wide fields also provide con- 
tinuous photometric coverage of tens of thousands of fainter 
stars, along with atmospheric transmission, cloud structure 
and auroral measurements. The initial survey demonstrated 
multi-color photometry on a;70,000 stars, showed that the 
images can be precisely astrometrically and photometrically 
calibrated, and demonstrated that the systems can achieve 
millimagnitude-level photometry. 

The paper is divided as follows. In section|2]we evaluate the 
expected transit detection efficiency at the Ellesmere Island 
site. In section[3]we describe the system design and construc- 
tion, including details on the thermal protection systems used 
to keep the optical elements clear of snow and ice, allowing 
continuous unattended operation in the Arctic conditions. In 
Section H] we describe the data collected during our test cam- 
paign, along with the data-reduction pipeline we developed 
for the systems. Section [5] details the image quality and pho- 
tometric performance of the cameras, and section|6]describes 
some initial astrophysical results from the dataset. Section |7] 
concludes with a discussion of concepts for much larger arctic 
surveys based on these camera prototypes. 

2. TRANSIT DETECTION EFFICIENCY AT ELLESMERE ISLAND 

Located on a 610-m ridge at the tip of the of the Fos- 
heim Peninsula, the PEARL facility - often referred to as 
the "Ridge Lab" - is accessed by a 15 km long road from 
Eureka, a research base operated by the meteorological di- 
vision of Environment Canada, the weather service of the 
Canadian government. The Eureka base is open to air traffic 
year-round, allowing astronomers to visit the site for observ- 
ing runs throughout the winter. 

Each year, the sun is continuously below the horizon at the 
Ridge Lab site for 3096 hours (129 days), and is continuously 
lower than -12° 1224hours (51days). Fall and spring "shoul- 
der seasons," when the sun is up for short periods, contribute a 
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Fig. 2. — A comparison of the probability of detecting a transiting planet 
in the field of a year-long survey, at a mid-latitude observatory and at an 
observatory at 80 degrees North. Because our exoplanet search is targeted 
at bright star's, we do not include the effects of moon brightness. The metric 
displayed here does not include the improved number of planet detections 
available if multiple fields around the sky are chosen for observation from 
the mid-latitude observatory; however, this cannot improve the numbers of 
long-period planets detected from mid-latitude sites by more than a factor of 
2-3 at best. To remove the effects of aliasing, each point in the curves is the 
average efficiency over a period range of one day. 

further 1,441 hours of dark time, for a total of 2,665 hours dur- 
ing which the sun is below -12° each year. The continuous 
darkness during the winter months allows continuous moni- 
toring of a field, eliminates day aliases in the measurement 
of variable objects, and improves the detection probability of 
a repeating event such as a planetary transit ( jPont & Bouchy| 
2005~) |Daban et al.|2010[|Crouzet et al.|20T0) i. 

To quantify this improvement, we performed Monte-Carlo 
simulations of the detection efficiency of a transiting planet 
survey as a function of orbital period for both a mid-latitude 
site and the Ridge Lab site (Figure 151. In the simulation, we 
assume that the telescopes only observe when the Sun is be- 
low an altitude of -12° and when the target is above an air- 
mass of 2.0. We remove nights in which the available obser- 
vation length is less than three hours, as is difficult to estab- 
lish a useful photometric baseline for transit detection in that 
time. We include weather-time-loss in each simulation, with 
24-hour periods of bad weather starting at randomly-chosen 
times (with overlapping bad-weather periods being allowed). 
We simulated correlated weather outages by assigning a 25% 
probability for a bad-weather period being immediately fol- 
lowed by another period of bad weather (and so on for subse- 
quent days). The weather simulation parameters were chosen 
to approximately reproduce typical weather outage patterns at 
mid-latitude sites (with a total of 35% of observing hours af- 
fected by weather), but provide a conservative estimation for 
the Arctic site, which often sees periods of hu ndreds of hours 
of clear skies in the winter) Steinbring et al. |2012[| . For the 
secure detection of an exoplanet candidate we require the de- 
tection of three transits. 

As shown in figure [2] Polar sites offer much improved de- 
tection efficiency for longer period planets, although we note 
that a mid-latitude site can improve the number of detected 
shorter-period planets by observing several fields throughout 
the year. Circumpolar surveys relying upon precision pho- 
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TABLE 1 

The specifications of the AWCam systems 
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Fig. 3. — A block diagram of the camera and data recording hardware. The 
cameras were placed in convenient locations on the Ridge Lab roof, with an 
approximately 100' cable run to the networking and data storage hardware in 
the PEARL control room. 

tometry or astrometry can also benefit from the near-zenith lo- 
cation of the celestial poles. At the Ridge Lab site an object at 
80 degrees declination can be continuously tracked through- 
out the winter months at an airmass range of 1.00 - 1.06, a 
small range of motion that can reduce atmosphere-induced 
systematic errors for many astronomical programs. 

3. INSTRUMENT DESIGN 

The AWCam prototype systems were required to operate 
autonomously in temperatures as low as -45°C, to avoid or 
remove ice and snow deposition on their optical surfaces, and 
to survive storm events including blowing snow with wind 
speeds as high as 40m/s. For these reasons, we designed the 
systems to be both simple and extremely robust. 

The cameras are based on a fast Digital Single Lens Reflex 
(DSLR) camera lens mounted to a 16MPix CCD camera. The 
wide-field-imaging system is aligned on the North Celestial 
Pole, and we take short 10s exposures to avoid the need for 
sky-rotation tracking. Continuous photometry of a circular 
area of several hundred square degrees around the Pole can 
be obtained. In the 2012 campaign we tested two systems, 
identical except in their lens and filter choices. The hardware 
and survey characteristics are summarized in Table [T] and an 
overview of the hardware is shown in Figure [3] 

The two camera systems were mounted on the observing 
platform of the Ridge Lab building, with Ethernet and AC- 
power connections to a warm control room within the build- 
ing. The control room also contained a network power switch 
for remote rebooting of the systems, a network switch, and a 
server for data storage. 

3.1. Camera assembly 

Each camera system is based on an f/1.2 Canon EF-mount 
DSLR lens, with focal lengths chosen to provide fields of view 
of hundreds of square degrees, and sufficiently small pixels 
(in the tens-of-arcseconds range) to allow precision photome- 
try of my < 10 stars without significant crowding effects. The 
lenses image onto Kodak KAF- 16803 CCDs packaged in a 
Finger Lake Instrumentation (FLI) PL- 16803 camera. A FLI 
CFW-4-5 filter wheel provides five filter positions for 50mm 
square filters. The Canon lens' built-in focus hardware can- 
not be easily accessed without connection to a Canon camera, 



Survey characteristics 



Pointing 

Survey dates 

Survey length (total) 

Survey length (dark and clear) 

Data collected 

CCD Hardware 



North Celestial Pole 

14 February 2012 - 21 February 2012 

152 hours 

98 hours 

44,583 images (1.36 TB) 



CCD 

Peak CCD Quantum Efficiency 
Pixel size 
Readout time 

85mm camera 



4096 2 front-illuminated (KAF- 16803) 

59% 

9/im 

4s 



Camera lens 

Field dimensions 

Continuous-coverage field 

Pixel scale 

Image quality 

Filters 

50mm camera 



Canon EF 85mm f/1.2L II USM 
25.4 X 25.4 degrees 
504 square degrees 
22.3"/pixel 

2-5 pixel FWHM over entire field 
Clear, g, r, i, z 



Camera lens 

Field dimensions 

Continuous-coverage field 

Pixel scale 

Image quality 

Filters 



Canon EC 50mm f/1.2L USM 
40.8 X 40.8 degrees 
1295 square degrees 
35.9"/pixel 

2-5 pixel FWHM over entire field 
Clear, g, r, i, z 



and manual focusing of these very fast lenses did not provide 
sufficient image quality. To allow precision focusing under 
remote computer control we added an FLI Atlas focuser to 
the optical system. The camera systems are controlled over 
a USB link by a compact Linux-based server located inside 
their thermal enclosures. The completed imaging system is 
compact, measuring a total of 230x230x175mm and weigh- 
ing 6.5kg. 

3.2. Enclosure, window and thermal system 

To ensure survivability of the electronics and other compo- 
nents, and to provide sufficient heating to remove snow and 
ice, the camera thermal enclosures were required to main- 
tain an internal temperature of 5C in temperatures lower than 
-40°C and in wind speeds up to lOm/s. We designed heated 
enclosures based on the Pelican 1610 case, an equipment 
shipping case which consists of a waterproof polypropylene 
copolymer enclosure. We drilled a 3.5-inch diameter hole on 
the upper surface for the camera aperture; a 3mm-thick AR- 
coated glass window protected the camera lenses. Two small 
holes equipped with bulkhead connectors provide inputs for 
AC power and Ethernet connections. The entire case internal 
surface was insulated by one-inch-think neoprene/vinyl/buna- 
N foam rubber. 

A high level of interior temperature stability was required to 
provide sufficient focus stability for the small depth of focus 
of our F/1.2 lenses. We based the case thermal control sys- 
tem on the heating and air circulation system of a Dotworkz 
Systems D2-RF-MVP CCTV enclosure which was used in an 
earlier camera prototype. We modified the control electronics 
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TABLE 2 
February 2012 observations 



Fig. 4. — One of the camera enclosures mounted at the Ridge Lab and 
aligned to the North Celestial Pole. Note the ice around the sides of the case 
and its surroundings, while the optical window is kept clear of accumulation 
by the heating system. 

to allow control by an Omega CN743 precision temperature 
controller equipped with an Omega HTTC36-K-14G-6 hol- 
low tube thermocouple probe. The heating system provides 
a total of 60W when fully activated, adding to the 70W pro- 
duced by the cameras, computers and other electronics in the 
enclosures. 

Pre-campaign testing of the systems was performed in a 
freezer at -27°C, with wind simulation provided by small 
fans. We found the enclosure's thermal control systems could 
maintain temperature inside the cases to ±0.5°C in conditions 
approximating those found at the Ridge Lab site. On-site, 
the systems were able to maintain the nominal 5°C internal 
temperature in all conditions encountered (-30+5°C), except 
when the wind speeds exceeded lOm/s. 

Precision photometry requires keeping the optical window 
of the systems free of snow and ice during arctic conditions. 
Our testing in the Arctic, including days with high winds, 
blowing snow and high levels of "diamond-dust" (fine ice par- 
ticles) demonstrated that the window was kept clear by the 
heating system in all conditions encountered (Figure ffl. The 
gentle heating allowed the snow and diamond dust to subli- 
mate off the camera enclosures rather than melting, reducing 
ice buildup around the instruments. 

3.3. Camera Mounting and Alignment 

The cameras were aligned on the North Celestial Pole by 
rotating the enclosures to face North and tilting them to the 
required angle. An initial azimuthal alignment was measured 
using a North-South line pre-determined via a handheld GPS 
receiver, and the system tilt was fixed using a digital incli- 
nometer. After initial pointing, fine alignment was achieved 
by monitoring the position of the North Celestial Pole in the 
acquired images. We used ratchet straps and shims to mount 
the enclosures onto the grillwork on the Ridge Lab roof, a de- 
sign which proved stable for all the encountered conditions. 

3.4. Camera Control Software & Data Storage 

The camera systems were controlled by custom software 
running on Ubuntu-Linux-based computers inside each en- 
closure. The software handled initial hardware tests, autofo- 
cusing, routine camera operations, and automatic transfer of 
acquired data to the data server. 

The autofocus routine was based on code used to fo cus the 
Palomar Transient Factory camera (Law et al. 2009). First, 
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the camera steps through a range of focus positions. Star-like 
objects in the field are ex tracted and have their F WHMs mea- 
sured using SExtractor (Bertin & Arnouts|1996| l. The result- 
ing list of objects is then filtered to exclude extended objects, 
cosmic rays and other non-point-sources, a median FWHM is 
derived for each focus position, and a parabola is fit to the re- 
sulting focus position / FWHM curve. The DSL R len ses pro- 
duce a changing PSF across the field (see section |5TT| ), includ- 
ing field curvature which can interfere with simple autofocus 
routines because different parts of the field can be in focus at 
different positions. We found that restricting the focus rou- 
tine to use only the central 80% of the image and rejecting 
stars with high (>1.5) ellipticity gave acceptable focusing for 
targeted zone of continuous coverage in the center of the field. 

The two camera systems each generated approximately 
8GB of image data per hour. We transferred the data im- 
mediately after acquisition onto two 2TB hard disks inside 
the Ridge Lab. The satellite Internet connection at the site 
was sufficient only for test and monitoring images to be trans- 
ferred during observations; the full survey dataset was stored 
on USB external hard disk drives for later analysis. 

4. OBSERVATIONS AND DATA REDUCTION 

Table [2] details the data taken during our February 2012 
campaign. Because the campaign was undertaken towards 
the onset of arctic sunrise, twilight conditions were encoun- 
tered for a fraction of the day, allowing a total of 16 hours 
of dark-sky observing in each 24-hour period. The cameras 
operated continually from 14 February 2012 to 21 February 
2012, a total of 152 hours of which 98 hours was dark. Thin- 
to-moderate clouds were experienced for approximately half 
the run, resulting in 0.5-2 magnitudes of extinction and in- 
creased scattering during those periods. 

To test different survey methods, we operated the two cam- 
eras systems in slightly different ways. The 85mm camera 
swapped sequentially between five filters between each expo- 
sure, providing near-simultaneous multi-colour photometry. 
The 50mm camera instead operated in a single fixed r-band 
filter. The r and i bands of the 85mm camera had poor im- 
age quality due wedge angles in the filter glass; we continued 
observations in those filters for sky-brightness measurements 
only. 

The cameras operated with an open-shutter efficiency of 
71% for the 50mm camera and 65% for the 85mm camera, 
which had slightly increased overheads because of its fil- 
ter changes. In total, the cameras collected 44,583 images, 
1.36 TB of data. 

4.1. Flat fielding with very-wide-field cameras 

Flat fielding with a fixed, autonomous camera system is po- 
tentially challenging. We tested four flat fielding methods to 
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Fig. 5. — Image cutouts from the 85mm-lens camera system, a 10s exposure in the g filter. Left: a central region from the image; Right, top-to-bottom: three 
small cutouts showing representative PSFs at the center, halfway to the edge, and the edge of the field respectively. 



inform subsequent instrument and survey design - conven- 
tional twilight flats, superflats generated from our full dataset, 
imaging a flat white screen held above the cameras, and using 
thick clouds as flatfields. 

4.1.1. Conventional twilight flats 

Flat fields taken during twilight are susceptible to large- 
scale gradients caused by multiple scattering. These ef- 
fects can generally be neglected for narrow-angle astronom- 
ical imaging, as long as extremely high precision is not re- 
quired. In contrast, our very-wide-field twilight flats showed 
unacceptably large gradients. This, coupled with the several 
months during which twilight is not available at circumpolar 
observing sites, makes twilight flats unsuitable for these cam- 
era systems. 

4.1.2. Superflats 

The major contributors to the sky brightness measured by 
superflats generated from data taken in dark time are airglow, 
zodiacal light, and starlighj^] Zodiacal light and airglow vary 
across the sky; at a zenith angle of 12.5 degrees (the edge 
of our narrowest field) the airglow is approximately 2.5% 
brighter than at the zenith, while the zodiacal light varies by 
approximately 10% acros s the range of ecliptic latitude cov- 
ered by our cameras (Benn & Ellison 1998). Because of our 
cameras' large pixels, the starlight contribution is larger than 

5 The Ridge Lab is located near the North Magnetic Pole, ensuring that 
most auroral emission is over the Southern horizon (Figure[TJ. 



found in narrow-field astronomical images, and the repeti- 
tive nature of the sky rotation complicates the generation of 
precision superflats. For these reasons, we found that it was 
extremely difficult to generate precision superflats from our 
dataset. 

4.1.3. Screen flats 

"Dome" flats generated from screens in front of the camera 
are not susceptible to the problems encountered by twilight 
and super flats, although care is required to ensure a uniform 
illumination pattern. Achieving this is simple for our fast, 
small lenses. We found that a white card held in front of our 
cameras during late twilight was sufficient to generate high 
quality flats (at a level sufficient for millimag photometry) for 
both our cameras and all filters. 

4.1.4. Cloud flats 

Without a dome on which to place the screen, screen flats 
cannot be easily taken by a standalone remote camera system. 
An automated system to place screens over the cameras would 
be complex and potentially unreliable. For these reasons, we 
also attempted to test the use of thick clouds as a substitute 
flat screen. The thickest clouds during which we obtained 
data, with an extinction of ~3.5 magnitudes, still had enough 
stars to complicate the generation of flat fields. Comparing 
the "cloud-flats" to the screen generated flat, we also found 
variations in the sky brightness at the 1 % level on scales of 5- 
10°, likely induced by structures within the cloud layer. This 
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could possibly be removed by averaging of a long-term cloud- 
covered dataset, but the cloudy data we collected was not of 
sufficient length to test this. 

4.1.5. Reliable Wide Field Camera Flats 

Twilight and cloud-based fiats appear to only be suitable 
for rough flat-fielding with these wide-field camera systems. 
Although the sky-gradient effects of scattering, airglow and 
zodiacal light could removed at some level using their pre- 
dicted across-field variations, it is difficult to generate clean 
flats with the stars in the field executing repetitive motions that 
puts them in the same pixels at the same time each day. Cloud- 
generated fiats showed percent-level spatial changes that will 
be difficult to accurately remove without a long period of 
cloudy weather, and there is a possibility that the brightness 
distribution of low-level clouds could be systematically biased 
by local topographic features. 

For these reasons, we elected to use the screen fiats for sub- 
sequent analysis. Remote operation of these systems is chal- 
lenging, although we note that the sealed optical system of 
our cameras and their very fast beams reduce the likelihood of 
changes in the fiats. We did not detect any fiat field changes 
over the course of our survey, at the one part in 1 ,000 level; it 
is likely that an occasional visit by a technician placing illu- 
minated screens in front of the cameras would be sufficient to 
achieve precision flat-fielding. 

4.2. Data reduction pipeline 

We implemented a Python-based data reduction pipeline to 
calibrate our images for photometry and astrometry. FITS 
files recorded by the cameras were first dark-subtracted and 
then flat-fielded using the methods described above . The 
sources in each fi eld were extracted using SExtractor ( |Bertin 
|& Arnouts||1996[ ) "mag-auto" apertures with modelling and 
subtraction of a background that was allowed to vary across 
the field. 

4.2.1. High-efficiency astrometric calibration of 
extremely-wide-field images 

The images produced by our lenses have a generally much 
higher degree of geometric distortion than those produced by 
standard astronomical telescopes. We found that automatic 
astrometric calibrations t ypically used f or wide-field s urveys 
such as astrometry.net (Lang et al. 2010) and SCAMP (Bertin 
2006) 1 led to unacceptably poor astrometric solutions (~100 
arcsecond median astrometric error compared to catalogs in 
some regions of the image), even with high-order polynomial 
corrections. We found that splitting the individual images 
into 10-30 rectangular subsections allowed SCAMP to find an 
acceptable astrometric solutio n within each subsection based 
on the UCAC-3 catalog (Zacharias et al. 2010 1. However, 
each subsection required tens of seconds to extract the sources 
and find an astrometric solution. This results in whole-image 
processing times of several minutes, unacceptable for images 
taken at least 4 times each minute in two cameras. 

To improve the processing times we developed a model to 
directly relate camera pixel positions and observation times 
to sky coordinates. Because the cameras' optical axes were 
slightly mis-aligned with respect to the Celestial Pole, we 
modelled lens distortion terms in polar coordinates around 
one "optical" axis, and then the conversion to RA and Dec- 
lination around a separate "celestial" axis. We modelled the 
lens distorti on in the f ollowing way, based on the models first 
described in Conrady| ( |1919| l and |Brown| (T966 ): 
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Fig. 6. — The photometric precision achieved for the 5,283 bright stars in 
the 50mm camera field of view. The red points show the RMS variations in 
the measured photometry across the 2801 individual photometric data points 
taken for each star over 12 hours of operations. The blue points show the 
stability when binned down to 10-minute (60 datapoint) chunks. 



D= J(x-x c ) 2 + (y-y c ) 2 



(1) 



Xlenscorr = (x~ X c )(l + aD 2 + bD 4 + cD 6 ) 

+2d(x - x c )(y - y c ) + e(D 2 + 2(x - x c ) 2 )) (2) 



yienscorr = C ~ JcXl + aD 2 + bD 4 + cD 6 ) 

+d(D 2 + 2(y- y c ) 2 + 2e(x - x c )(y - y c )) (3) 

where x,y are the measured stars positions in pixels, x c ,y c 
are the camera optical axis pixel coordinates, and the terms a- 
e are fit to the measured astrometry. To measure the model co- 
efficients and optical axis coordinates we took 20 frames dis- 
tributed throughout the nights, fully calibrated the frames us- 
ing astrometry.net and SCAMP, and used a downhill simplex 
minimization to fit the coefficients to the calibrated astrome- 
try. After correction of the pixel coordinates, the sky coordi- 
nates of the stars are calculated by converting the de-distorted 
pixel positions to sky coordinates with a tangent plane projec- 
tion and then adding the sky rotation term corresponding to 
the observation epoch. 

After initial modelling, we found that a fixed lens model did 
not adequately fit the data. An additional shift was required 
to model a slow change in the angle between the optical and 
celestial axes. The amplitude of the modelled shift is consis- 
tent with a 100/im height change in one edge of the camera 
support structure, and was maximal around midnight. We at- 
tribute this to the effects of changing temperature during the 
semi-twilight nights. 

Verifying the full astrometric solution against a second set 
of calibrated images, we found that just based on the time of 
observation and the pixel position of the stars, sky positions 
were calculated to «20 arcsec (<1 pixel) accuracy across the 
entire image, with improved performance in the image cen- 
ters. The full astrometric processing for the ~ 10,000 stars in 
each image takes a fraction of a second. 

4.2.2. Efficient millimagnitude photometry with 
extremely-wide-field images 
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To obtain photometric zeropoints we modified the differen- 
tial aperture photometry pipeline originally developed for the 



Palomar Transient F actory ( Law et aLp 009 ) M-dwarf plane- 
tary transit search ( [Law et al.|2011| l. We associate an initial 
"rough-guess" photometric zeropoint with each image, and 
using them generate light curves for each source. We then 
vary the zeropoints, using a simulated-annealing algorithm to 
minimize the averaged noise-weighted RMS variability in all 
the light curves. The algorithm converges on the set of ze- 
ropoints which provide the most stable photometry for the 
largest number of sources in the field, and is thus robust to 
bias from astrophysical or instrumental variability in individ- 
ual stars. 

The slight mis-alignment betw een the cameras' optical axes 
and the Celestial Pole (Section |4.2.1[ ), and the small focal 
plane tilt (Section |5.1| l, leads to a slow change in size and 
shape of the stars' PSFs as the sky rotates. We found that this 
was enough to induce few-percent-level changes in the post- 
zeropoint-calibration stellar fluxes, for stars at the edges of the 
field where the PSF changes were largest. The resulting sys- 
tematic flux error affects each star in a particular region of the 
image in the same manner, and so can be removed with stan- 
dard methods designed to rem ove light curve systematics. We 
use 15 iterations of SysRem (Tamuz et al. 2005 1 to model and 
remove this systematic error for the ensemble of light curves. 

After photometric calibration and systematics removal, 
we consistently achieved 1 -percent-level photometric perfor- 
mance in individual images over 12-hour timescales, and few- 
millimagnitude-level photometry when binning datapoints 
from 60 individual exposures together (Figure |6j. 

The expected single-exposure scintillation noise for at- 
zenith observations is approximately 6 and 8 millimagnitudes 
for the 85mm and 50mm lenses respectively ( |Dravins et al.| 
1 1 99 8 1 |Young| 1 967) . The scintillation noise estimates depend 
on the particular observing conditions at the site (to within 
~50%; |Young| [i967 ), and we thus conclude that the cam- 
eras are achieving very close to scintillation-limited photo- 
metric precisions. We evaluate the photometric performance 
for planet-searc h pu rposes, including the level of correlated 
noise, in section [53| 



5. PERFORMANCE 

In the below subsections we detail the performance of the 
camera systems in the context of the requirements of a bright- 
star transiting exoplanet survey. 

5.1. Image Quality 

The lens point spread functions (PSFs; Figure Bl) are well 
matched to the camera's pixel size in the center of the field, 
with round shapes and full width at half maxima (FWHMs) 
of approximately two pixels. In both cameras, the PSFs in- 
creased in width towards the edges of the field (Figure [71, 
although at a slow enough rate to provide useful images up 
to the edges of the zone in which stars are continuously cov- 
ered. A small tilt in the focal plane (~5yum across the CCD) is 
visible in each camera's FWHM map, produced by a combi- 
nation of the camera mounting hardware and the tolerance in 
the wedge angle of the filter glasses. The tilt did not signifi- 
cantly affect the image quality within the continuous coverage 
zone. 

5.2. Camera Sensitivity 

The cameras achieved their design sensitivity, with a typi- 
cal lOs-exposure 5-sigma point-source limiting magnitude of 
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Fig. 7. — Maps of the FWHM of stars recorded in the standard survey bands 
for the two cameras: the g filter using the 85mm F/1.2 lens (top) and the r 
filter with the 50mm F/1.2 lens (bottom). The circle marks the zone for which 
the camera PSFs were optimized, where stars have continuous photometric 
coverage as the sky rotates. 

m g = 12.6 for the 85mm lens. We measured the limiting mag- 
nitudes using s tars observed in a Sloan Digital Sky Survey 
SEGUE field ( |Yanny et al.|2009b which was included in both 
our cameras' fields of view and allows us to approximately 
match filter responses for photometric calibration. At its high- 
est point, the SEGUE field covers declinations from 82.3° to 
85.1°. Our cameras have somewhat lower image quality and 
increased vignetting at those altitudes than at the center of the 
fields. Correcting for these effects, the true limiting magni- 
tude at the field centres is improved to m g « 13.6. 

Figure [8] shows the relative vignetting of our cameras as a 
function of field position. The 85mm lens reduces the light 
collection to approximately 50% of its maximum value by the 
edge of the camera's zone of continuous photometric cover- 
age, while the 50mm lens has a more pronounced vignetting 
giving ^35% sensitivity at the edges of the chip. 

5.3. Correlated noise performance 

As shown in Figure [6] binning our single-exposure pho- 
tometry allows few-millimagnitude precision for thousands of 
stars in the our fields of view. However, correlated ("red") 
noise can greatly reduce the efficiency of planet detection 



8 



N.M. Law et al. 





Fig. 8. — The vignetting of the fields of the 85mm (left) and 50mm (right) 
lenses: the Hatfield images from the cameras. A sensitivity of at least 50% of 
the maximum value is maintained throughout the continuous-coverage zone. 
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Fig. 9. — The photometric stability (RMS compared to median magnitude) 
achieved by the 50mm camera (top) in the r filter and the 85mm (bottom) 
camera in the g filter on a typical my = 7.0 star, which is bright enough for 
its photometry to be limited by scintillation and systematics. We binned in- 
dividual lOs-exposure points to achieve 2-3 millimag stability over effective 
exposure times of 5-10 minutes. The expected behaviour as exposure time 
increases, with the photometric precision improving with the square root of 
the total exposure time, is a close match to the measured performance. 



Smith 



even with apparently good photometric stability (e.g. 
|et al.|2006| l. 

We evaluated the correlated noise performance of our cam- 
era systems by searching for unexpected behaviour when bin- 
ning the photometric data points to different timescales; a sys- 



tem with a high degree of correlated noise would show a pre- 
cision improvement with binning which is slower than the ex- 
pected sjNdata points- To remove the possibility of obtaining 
spuriously good performance estimates by over-fitting during 
the systematic noise reduction steps, we used a subset of our 
full data reduction pipeline, without the final photometry and 
systemic noise removal steps. We repeated the below tests 
with our full pipeline and obtained identical results. 

Using the SExtractor extraction of the sources in our 
astrometrically-calibrated images, we performed differential 
aperture photometry for four isolated stars near the center of 
our fields. We selected the target and three calibration stars to 
have magnitudes in our exoplanet survey range (m v = 7- 9), 
but with care to keep their peak fluxes below ~20, 000 ADUs 
to avoid non-linearity due to our CCD's anti-blooming system 
(which engages at signals of -^32,000 ADUs). Within each 
night, we fit a second-order polynomial to the light curves 
to remove long-timescale residual systematic variations (typ- 
ically at the 5-millimagnitude level). The final photometric 
stability was measured from the RMS variability in the result- 
ing light curve. 

Similarly to the full pipeline, we achieved an RMS preci- 
sion of 1.0% in individual 10s exposures on the 85mm cam- 
era, and 2.4% in r-band exposures with the 50mm camera. 
The ratio of photometric performance for the two lenses is 
roughly as expected for measurements limited purely by pho- 
ton and scintillation noises, suggesting only low-level extra 
sources of photometric error. 

Binning the photometry to 10-min exposures improved 
the photometric precision to 2-3 millimagnitudes (Figure [9]), 
compared to an expected scintillation-limited precision of «2 
millimagnitudes. In both cameras, the photometric precision 
improves the expected rate of sjN 'data points until at least 100 
frames are binned together. The length of our dataset from 
this campaign precludes a meaningful analysis of the effects 
of correlated noise on longer timescales, but the close adher- 
ence to the ^Ndata points relation is encouraging. 

6. DETECTED VARIABLE SOURCES 

We performed a search for variable stars and other astro- 
physically varying sources, both visually searching for obvi- 
ous periodic sources, and testing automatic variable-source 
detection in a blind search for stellar variability and transit 
signals. 

6.1. Bright periodic sources 

Because our targets are bright (my <~ 10), essentially all 
the high-amplitude variables in the field are already known. 
In an initial visual characterization of each light curve with 
evidence of a large amount of variability, we re-detected the 
Delta-Scuti star V377 Cep, and detected eclipses for the beta 
Lyr eclipsing binaries EG Cep and AZ Cam, the W UMa 
eclipsing binaries FN Cam and RZ UMi, and the Algol-type 
Binaries W UMi, TY UMi, AY Cam and SV Cam. In Figure 
10 we show a two-camera, three-colour, 15-hour light-curve 
or the primary eclips e of W UMi, a my = 8.5 Algol-type 
eclipsing binary (e.g. Sahade|1945 1. Three primary and sec- 
ondary eclipses of this 1.71 -day-period system were captured 
during our week-long observations. 

6.2. Stetson-J blind detection of variable stars 

We tested the ability to perform an automated stellar vari- 
ability and transit search with the cameras, using known 
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Fig. 10. — The 18 February 2012 primary eclipse of W UMi, followed over 
a period of 15 hours by both camera systems. The small grey points are 
individual 10s exposures from the 50mm camera in the r filter; the lines and 
points with error bars are binned photometry from each camera, with error 
bars estimated from the RMS scatter of points within the bin. The red points 
are near-simultaneous binned photometry taken by the 85mm camera in the 
z filter, binned over 2x more points because of the decreased sensitivity of 
the camera in the z filter. The increase in photometric scatter towards later 
times is due to the beginning of twilight; the z filter is more affected by this 
and cuts out earlier. The photometric precision of the binned points is 2-5 
millimagnitudes, depending on the filter and camera. 
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Fig. 11. — Stetson-J variability index for unbinned light curves of the 5,283 
stars in the 50mm-camera exoplanet-search sample. The Stetson-J index for 
the binned data is similar but with a different normalization. 



bright variable objects in our field to evaluate the detection 
efficiency. The Gen eral Catalogue of Variable Stars (GCVS; 
Samus et al. 2009) lists 119 variable stars within our field 
of view and magnitude limits (5.0 < my < 10.0). Most of 
these objects are rarely outbursting stars or long-period vari- 
ables which are unlikely to have produced detectable variabil- 
ity during our observations. However, 28 stars have <20 day 
photometric periods (and so a reasonable chance of undergo- 
ing measurable variability within our observations), and have 
light curves which were flagged as clean of instrumental ef- 
fects by our pipeline. All the known variables in the list have 
amplitudes which can be easily detected by our cameras. 

We based the v ariability search on the Stetson-J statis- 
tic (Stetson 1996}, which compares pairs of observations 
to search for variability which is correlated between data- 



points. Compared to a simple RMS -variability selection we 
found that Stetson-J is much less sensitive to the main source 
of false-positive variability in our dataset: few-percent-level 
blending between closely separated sources, which affects ap- 
proximately 2% of the stars in the field. Most of the stars in 
the field move by a large fraction of a pixel each frame, and 
thus the blending between closely-separated stars can vary 
by large amounts (> 10%) frame-to-frame. The resulting 
mostly uncorrelated blending from frame-to-frame allows the 
Stetson-J statistic to exclude those sources as instrumentally 
rather than astrophysically variable. 

To detect a source as variable we set a minimum Stetson- 
J statistic of 0.14, a value selected to exclude 98% of the 
sources in our field (Figure [TT}. Seven (25%) of the GCVS 
variable stars were detectedas variable under this crite- 
rion. The light curves of the remaining GCVS objects were 
clean of obvious instrumental effects, but did not contain 
eclipses or other obvious astrophysical variability, suggesting 
the pipeline is efficiently and automatically detecting the ob- 
jects that underwent eclipses or other variability during the 
observation period evaluated. 

A visual inspection of the light curves of the other 1 19 high- 
variability sources revealed eight objects with signs of pul- 
sations or other sinusoidal variability at the <5% level; five 
objects with possible single eclipses with 5-20% depths (sug- 
gesting long-period systems); 43 sources with blending ef- 
fects from nearby stars; and 63 objects with long-timescale 
variability that we cannot usefully classify as instrumental 
or astrophysical in the current dataset. The planned full- 
winter Arctic datasets will allow confirmation of these pos- 
sible longer-period variable stars and enable a full evaluation 
of stellar variability in this field on timescales up to several 
months. 

6.3. Transit search 

We searched for individual planetary transits using the 
binned high-precision light curves shown in Figure[6] We cal- 
culated the expected Stetson-J transit signal by injecting sim- 
ulated transit events (1%, 2-hour-duration dips with ingress 
and egress periods) into our binned light curves. We found 
that individual transit events are securely detectable using 
the Stetson-J statistic for the brightest 1/3 of the stars in 
the dataset (my ^8.5). Typical 1%-depth transits produce 
Stetson-J values of 0.3 in the binned light curves of these 
bright objects, excluding the 97.5% of objects in the field 
which have lower values. A visual search of the 38 light 
curves with Stetson-J > 0.3 in the binned light curves revealed 
no likely transit-like events. The statistically-significant de- 
tection of transit signals around the fainter stars will be fea- 
sible with a longer dataset which allows the use of more 
sophisticated phased matched-filter-type transit-search algo- 
rithms that o btain mu ch improv ed detection significances (e.g. 
Kovacs et al.|2002"l|Tingley|2003| l. 



7. DISCUSSION & CONCLUSIONS 

We have demonstrated the operation of two wide-field cam- 
era systems in the High Canadian Arctic, achieving the image 
quality and photometric precision necessary for the detection 
of transiting exoplanets around bright stars. Each camera's 
dataset contains light curves for «70,000 stars, with «5,000 
in each camera having sufficient SNR and being sufficiently 
uncrowded to search for transiting exoplanets. 

In a dataset from a full arctic winter, the probability of de- 
tecting at least three transits of planets in month-long orbits 
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around solar-type stars is «70% (Figure[2]i, while the geomet- 
ric transit probability for those planets is a few percent. With 
few-millimag photometric precision the systems could detect 
exoplanets as small as Neptune transiting solar-type stars. The 
Kepler-derived few-percent occurrence frequenci es for plan- 
ets in that size and orbital period range ( Howard et al.pOlF 



suggest that the AWCam systems, when operated during an 
entire arctic winter, will be capable of detecting up to five 
transiting exoplanets around stars brighter than mv=9.5. 

7.1. A concept for a 10,300 square degree Arctic camera 

Although wide field, the cameras described here only cover 
3-6% of the sky continuously accessible from the Ridge Lab 
in winter months. Larger sky areas can be covered with mul- 
tiple cameras, but this requires extremely short exposures or 
a tracking mount to avoid star trails. One such concept un- 
der dev elopment, the Co mpound Arctic Telescope System 
(CATS, 



Law et al. |2012 1, places 19 cameras similar to the 
ones described in this paper in a hemisphere which rotates 
once per day. This arrangement would cover a total of 10,300 
square degrees from one rotating mount, providing full cov- 
erage of the Polar sky down to 30° degrees declination. The 
system is designed to be capable of monitoring 90,000 stars 
brighter than 10 th magnitude, with millimagnitude precision 
in each few-minute exposure. A further ~10 6 stars can be 
covered with millimagnitude precision with a 30-minute ca- 
dence, and stars brighter than my ~ 17 would be monitored 



with SNR>5 every few minutes. The camera has only one ma- 
jor moving part to provide the mount rotation, and so can be 
constructed to be sufficiently reliable for autonomous opera- 
tions in the Arctic. Coupled to a small spectrograph-equipped 
telesco pe, such as an upgraded Dunlap Institute Arctic Tele- 
scope (Law et al. 201 1| [2012), the system would provide an 
integrated, rapid High Arctic survey and follow-up facility. 
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